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Abstract 

The recently proposed probability representation of quantum mechanics is generalized 



' to quantum field theory. We introduce a probability distribution functional for field con- 



figurations and find an evolution equation for such a distribution. The connection to the 
time-dependent generating functional of Green's functions is elucidated and the classical 
limit is discussed. 

The use of statistical methods for describing quantum physics gives the opportunity to 
describe classical and quantum phenomena in a unified approach. Since the beginning of 
^ - quantum mechanics there have been attempts to understand its nature in a classical-like con- 

■ text, namely to describe quantum states in terms of a classical distribution of probability. 

It is this philosophy which inspired the so-called quasi-probability distribution functions of 
Wigner, Husimi, Glauber and Sudarshan The original goal was not completely achieved 
(the above distribution functions are not always positive defined or they do not describe mea- 
surable variables) until Cahill and Glauber in Q introduced a class of distribution functions, 
known as marginal distribution functions (MDF), which enjoyed all the properties of a density 
of probability. Nevertheless, it was realized only recently Q that quantum mechanics could be 
described entirely in terms of a distribution of such a family, suitably defined for a random vari- 
able, which we will specify below. In Q a consistent scheme has been proposed, the so-called 
probability representation, which has been shown to be completely equivalent to the ordinary 
formulation. Quantum states are described by a distribution of probability, the MDF, and the 
time evolution by an integro-differential equation for the MDF. Invertible relations have been 
established between the MDF and the density matrix |^, ^ and between the Green's functions 
of the related evolution equations [||. 
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In this framework classical and quantum phenomena, both statistically described, only 
differ by the evolution equations of the distributions of probabilities for the relevant observ- 
ables. The quantum evolution equation, of Fokker-Plank type, is seen to reduce to Boltzmann 
equation for the classical distribution of probability when the classical limit is considered. 

In this letter we generalize the probability representation first to the case of interacting 
particles, then to non relativistic quantum field theory. We consider a system of interacting 
oscillators which describe, when the limit to the continuum is performed, a self-interacting 
scalar field theory with generic self-interaction potential. We introduce the notion of MDF for 
the quantum state of such system and derive the evolution equation both for the discrete and 
continuous cases. The MDF is seen to be a distribution of probability with the same arguments 
used for 1-d quantum mechanics. 

Interestingly similar ideas have been developed by Wetterich in in connection with the 
approach to equilibrium in non-equilibrium quantum field theories. There an evolution equa- 
tion is found for a suitably defined time-dependent generating functional 0. We establish the 
connection between our MDF and Wetterich's generating functional in terms of a (functional-) 
Fourier transform. 

In section 1 we briefly review the probability representation of quantum mechanics for 
the simple case of a one-dimensional quadratic Hamiltonian. The derivation of the evolution 
equation for the MDF is performed in detail in a slightly different manner from its original 
derivation but more suitable for generalizations. In section 2 we consider a quadratic 
Hamiltonian describing N interacting particles. We define the MDF and show that this is a 
well defined probability distribution. We then find the evolution equation. In section 3 we 
consider a scalar field theory with self-interacting potential, which may be seen as the limit 
to the continuum of the previous model. We define the MDF which is now a probability 
distribution functional, and derive its exact evolution equation. In section 4 we derive the 
evolution equations for the above mentioned systems, directly in terms of the Fourier transform 
of the MDF, the quantum characteristic function. We show then how our results may be 
connected to those found in p. 



1 The Probability Representation of Quantum Mechanics 

The MDF of a random variable X was introduced in as the Fourier transform of the quantum 
characteristic function x(^) =< ^^'^^ >i to be 

wiX,t) = — I dke-'^^ < e'^^ > , (1.1) 
27r J 

where X is the operator associated to X, < A >= Tr:{pA), and p is the time-dependent density 
operator. It is shown in |2| that w{X, t) is positive and normalized to unity, provided X is 
an observable. This theorem may be easily proven taking for simplicity p to be the density 



operator for a pure state. Then, evaluating the trace in (1.1) on eigenstates of the operator 
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X, it can be verified that ( |1.1D yields w{X,t) 
unity. 



p{X, X, t) which is positive and normahzed to 



We recaU that the quantum characteristic function is, up to factors of i, the generating 
function of the momenta of any order, for the probabiHty distribution of the operator X. 
Hence it plays in quantum statistical mechanics the same role as the generating functionals 
for the Green's functions in quantum field theory. In ref. |^ X is taken to be a variable of the 
form 

X = fiq + iyp, (1.2) 

where n, v are real parameters labelling different reference frames in the phase space, [i is 
dimensionless, while \y\ = {m~^\\t\. Thus, X represents the position coordinate taking values 
in an ensemble of reference frames. For such a choice of X it was shown that there exists 
an invertible relation among the MDF and the density matrix, respectively in Q for the 1-d 
case, and in Q| for the 2-d case. This relation was originally understood through the Wigner 
function: the MDF was expressed in terms of the Wigner function which is in turn related 
to the density matrix and viceversa. The evolution equation of the MDF was then found 
starting from an evolution equation for the Wigner function established by Moyal in This 
intermediate step in terms of the Wigner function is not necessary. We can directly invert ( |1.1D , 
when the variable X and the associated operator are given by (|l3). The evolut ion equation 
is then obtained (in the Schrodinger representation) by means of the Liouville equation for 
the density operator, in coordinate representation. In view of the subsequent generalization to 
degrees of freedom and to field theory, let us derive these results in some detail for a one 
dimensional system. Equation (^]^) is explicitly written as 



w{X,ij,,u,t) = — 
zvr 



dk 



— / dk 

2tt 



dZ p{Z^ Z — kuK)e 



-ik[X-/i{Z-kuh/2)] 



(1.3) 



The MDF so defined is normalized with respect to the X variable: j dXw{X, p,,u,t) = 1. 
Performing the change of variables Z' = Z, Z" = Z — kvTi we may reexpress the MDF in the 
more convenient form: 



w{X,p,u,t) 



1 



2'it\u\Ti 



p{Z',Z",t)exp 



.Z'-Z" 



X 



Z + Z' 



which can be inverted to 



^ ^1 

p{X,X',t) = \a\ I w{Y,p, — — — )exp 



ia [Y — p 



X + X' 



dZ' dZ" 



dpdV 



:i.4) 



:i.5) 



where a is a parameter with dimension of an inverse length. The density matrix is independent 
of a. In facts, using the homogeneity of the MDF, w{aX,ap,ai') = \a\~^w{X, p,^), which is 
evident from the definition, (|1.5| ) may be written as 



p{X, X', t) = j w{Y, p,X - X') exp 



i ( X + X' 



dp dY 



:i.6) 
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where the variables Y, ^ have been rescaled by a. It is important to note that, for ( |1.4D to be 
invertible, it is necessary that X be a coordinate variable taking values in an ensemble of phase 
spaces; in other words, the specific choices = 1, i/ = or any other fixing of the parameters jjl 
and V would not allow to reconstruct the density matrix. Hence, the MDF contains the same 
amount of information on a quantum state as the density matrix, only if Eq. (1.2) is assumed. 

We now address the problem of finding the evolution equation for the MDF, for Hamilto- 
nians of the form 



Using the Liouville equation 



(1.7) 



:ii 



and substituting into Eq. (1.4), we have 



2'n\v\'h 
X exp 



2m \ dZ^ 9Z'2 



Z-Z' 
vfi 



X - n 



Z + Z' 



{V{Z)-V{Z')) 
dZdZ' . 



p{Z,Z\t) 



:ii 



Integrating by parts and assuming the density matrix to be zero at infinity, we finally have 



w{X,fj,,v, t) 



1 d i 
m ov n 



V 



V 



^dX' dn 2 dX 



d ]^ d ivfi d 
^dx' dJl~~dX 
iuh d 



w{X,n,u, t) 



where the operator (^) ^ is so defined 



i^r'J f{Z)e^^'^^dZ = J l^e^i^^^dZ . 



(1.10) 



(1.11) 



This equation, which plays the role of the Schrodinger equation in the alternative scheme just 
outlined, has been studied and solved for some quantum mechanical systems [10|,[11|. The 
classical limit of (1.10) is easily seen to be 



w{X,fi,u,t) = \^^^ + uV' 
m ov 



\dX J dfi / dX 



wiX,^,^, t) 



(1.12) 



where V is the derivative of the potential with respect to the argument. Equation ( |1.12| ) may 
be checked to be equivalent to Boltzmann equation for a classical distribution of probability 
f{Q,P,t) , 

(1.13) 



dt m dq dq dp ' 



after performing the change of variables 

w{X,^Ji,v,t) = ^J f{q,p,t)e"'^''-^^''-'P^dk dq dp 



(1.14) 



Hence, the classical and quantum evolution equations only differ by terms of higher order in 
h. Moreover, for potentials quadratic in q, higher order terms cancel out and the quantum 
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evolution equation coincides with the classical one. This leads to the remarkable result that 
there is no difference between the evolution of the distributions of probability for quantum and 
classical observables, when the system is described by a Hamiltonian quadratic in positions 
and momenta. For this kind of systems, the propagator is the same [^j. Of course, what makes 
the difference is the initial condition. 



2 Generalization to degrees of freedom 

We consider now a system of N interacting particles sitting on the sites of a lattice (we choose 
it to be one-dimensional for simplicity). The Hamiltonian of the system is 

N ^2 

'^ = i:|^ + vm (2.15) 

We assume the masses to be equal to unity. We take the potential to be of the form: 

^ / 1 \ 
^(^) = E (^('^"^+1 - '^^)' + ^('1^^) (2.16) 

where a is the lattice spacing and U{qi) is the part of the potential which depends only on 
the position of the i-th particle. This specification is not essential for the purposes of this 
section, but it will become necessary for understanding the limit to the continuum, which will 
be considered in next section. The quantum characteristic function for the dimensional 
system may be defined as 

x(fci,...A:^)=<e^S.^»^» > , (2.17) 

where < A >= Ti^pA), and p is the density operator of the system. (In case there is no 
interaction between different sites of the lattice the density operator may be factorized and the 
characteristic function is just the product xi^i, ■■■kN) = ^iLiX{ki) •) 

Performing the Fourier transform of ( p. 17] ) the MDF is then given by 

w{X^,^,^,iJ^,t) = j^J dA:i...dA;jve^'S,'=-^» < e^E,'^^^^ > ; (2.18) 

where cr is a collective index. It may be shown that this is a probability distribution, namely 
that it is positive definite and normalized, provided Xi are observables. The proof goes along 
with the one-dimensional case. We first suppose that there is no interaction between different 
sites at some initial time to and we assume for simplicity that the system be in a pure state 
\tp >= \tpi > (gi... (gi \ipj\f >. Using the factorization property of the quantum characteristic 



function Eq. (|2.18D may be seen to reduce to the product w{Xcr,Pa, i^a,to) = I\i Pi{Xi, Xi,tQ) = 
p{X,X,t()), which is positive and normalized. Then Liouville equation guarantees that this 
result stays valid when the interaction is switched on. In analogy with the one-dimensional 
case we now introduce the variables 

Xi = piQi + i^iPi (2-19) 
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with Xi accordingly defined. 

Introducing the notation \Zu >= \Zi > > we rewrite ( |2.18| ) as 

= ^ / Ul,dhdZ,p{Z^,Z^ - fc^^^;,)e-'E;+.^.[^.-/^.(^.-'=.-.V2)]^ 
where p{Z„, Z'^) =< Zfj\p\Z'^ >. Performing the change of variables 



(2.20) 



Z\ — Zi , Z\ — Zi — hiVifi 



(2.21) 



we have 



X exp 



N 

n 

N ryl ryll 



Xj - pj- 



Z' + Z'' 



(2.22) 



which can be inverted to 



(2vr) 



N 



N 



Yl dfii dYi w{Y„,pa,X^ - X'^) exp 



i=l 



N 



(2.23) 

Once again, we recall that the inversion of (|2.22| ) is made possible by choosing the variables 
Xi as in (|2lg| ). We now use the Liouville equation (|l.8|) to get 



w{X„,p„,v„,t) 



N 



i=l 



1 



+ {V{Z,)-V{Z'^)) 



Y[[—dZidZl 



N 



Lj=l 



2m\dZ] dZf ^ 



(2.24) 



TV 



Integrating by parts and assuming the density matrix to be zero at infinity, we finally have 



N 



w{X^,fi^, v„,t) 



'dX„' dp^ 2 dX, 



dX„' dpa 2 dX„^ 

w{X„,^i„,v„,t) , (2.25) 



where the inverse derivative is defined as in ( 1.11 ). We report for future convenience the term 
containing the potential when explicitating the interaction between neighbours: 



'dX^' dfia 2 dXa 



V 



d , 1 d iU(jTi d 



in X - 

1=1 



dX„' dn^ 2 dX^ 
d f d d 



'dX„' dfi^ 2 dX„ 



'dX^' dp^ 2 dXo 
2—+ ' ^ 



dXi \dXi+iJ dm+i dpi dXi \dXi^i) dpi-i 



(2.26) 
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When considering the classical limit we have 



d 1 



d f d 



d d 

2 \ 

d^i dXi 



d 



1 d 



d 



+Vi 



(-'at'-' 



dXi.i 

d 



i-l 



d 



dfiaJ dX„ 



w{X^,li„,u„,t) 



(2.27) 



where Ui is the derivative of the self-interaction potential with respect to the i — th vari- 
able. Equation ( 2.27] ) may be seen to be equivalent to the Boltzmann equation as in the 
one-dimensional case. Moreover, Hamiltonians which are quadratic in positions and momenta 
yield the same evolution equations for classical and quantum probability distributions. 



3 Generalization to Field Theory 

We now consider a scalar quantum field theory described by the Hamiltonian 



H 



(fx 



6=1 



(3.28) 



d is the spatial dimension, while U{(j){x)) is the self-interacting potential, polynomial in the 
field (p. The Hamiltonian ( p. 28 ) is easily seen to be obtained by the discrete Hamiltonian ( p. 15 ) 
by taking the limit to the continuum (a — > 0) with the following rules: 



d'^x 



„-(d/2+l) 



7r(x) 
(gi+i - qi 



d(\){x) 
dxh 



(3.29) 



In analogy with the discrete case we introduce the field 

d(x) = ^i{x)^{x) + v{x)tt{x) (3.30) 

where n{x) = lima^o '^^^/^i) and v{x) = \mia^Qa~'^/'^Vi . 

The quantum characteristic functional, which now will play the role of generating functional 
for correlation functions of the fields, may be defined as 

X{k{x)) =< e^/'^'^ >= Ti (p(t)e*/'^'^ fcW-IW^ ^ (3 3^^ 

The functional Fourier transform of xi^i^))-, what we will call the marginal distribution func- 
tional (MD.F), still defines a probability distribution. This can be understood by recognizing 
that it is the limit of the MDF for the discrete A^-dimensional system considered in the previous 
section: 



w{Mx),^i{x),v{x),t) = I Pfce-^/'^(^)*(^)'^^x(A;) = limw;(X,,/x,,i/,,t) 



(3.32) 
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where Hi (2tt)^ ~^ ^ density matrix functional may be defined as the hmit of 

Eq. (|2]2|) to be 



dy 



. (3.33) 



Then, the evolution equation for the probability distribution functional is easily obtained by 
taking the Hmit of ( ^^ ): 

-1 X 



w{^{x),n{x),v{x),t) 



(fx 



5v{x) 5^{x) 



5(t{x) J 5fi{x) 



+ 



U 



U 



/ -(5 6 iy{x)n 6 

\6^{x)) 5fi{x) 2 6^{x) 



M{x)J dfi{x) 2 6^{x) 
The inverse functional derivative (;5^^) is so defined: 



w{^{x),fi{x),u{x),t) 



(3.34) 



<5$(x) 



Vk e 



i J k{y)'i>(y)dy 



Vk 



k{x) 



-i J k(y)^(y)dy 



(3.35) 



while the notation A[/(x)] stands for f{x + Ax) — f{x). Performing an expansion in powers 
of % the classical limit may be obtained as in the previous sections. 



4 The Quantum Characteristic Function as a Generating Func- 
tion 

In this section we discuss the connection between the probability representation described 
above both for quantum mechanics and quantum field theory and a slightly different point 
of view developed in Q, where evolution equations are found for a suitably defined euclidean 
partition function. There are two main ingredients in our approach: one is is the probabilistic 
interpretation for the distribution describing the observables, the other is the equivalence 
between the description based on the MDF and the conventional description based on the 
density matrix. The first aspect is guaranteed by the Glauber theorem which states that 
the Fourier transform of the quantum characteristic function associated to observables is a 
probability distribution. The second aspect, namely the invertibility of the MDF in terms 
of the density matrix, is achieved by introducing configuration space variables which take 
value in an ensemble of reference frames in phase space, each labelled by the two parameters, 
/i, V. Thus, the evolution equations which we have found ( (ITlOl) , (|]2|), (|]3|)), together with 



suitable initial conditions, completely characterize the state of the given quantum system. 
These equations assume a simpler form when their Fourier transform is performed. We have 



X{k, ^i, y,t)= i dXe^'^'wiX, fi, u, t) , (4.36) 
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with obvious generahzations to the N dimensional case and to field theory. For the one- 
dimensional quantum systems considered in section 1, the Fourier transform of Eq.( [L.ld ) yields 
an evolution equation for the quantum characteristic function itself: 



. , I 1 5 i 

' m ov n 



ik dji 



^ 1 d kuh 
Jkdn 2 



X{k,fx,u,t). (4.37) 



For the A^-dimensional quantum systems considered in section 2, the Fourier transform of 



Eq.(|2^) yields 



X{ka,Ha,l'a,t) 



1 d i 



V 



V 



19^ Ufjkcrfl 



_1 d_ 

ikfj d/j^a- 2 



Jka dfj,cr 2 

while Fourier transforming Eq.( 3.34| ) we have, for quantum field theory. 



(4.38) 



X{k{x),fj,{x),i^{x),t) 



cfix 



S 



m ov[x) 



2ihk{x)i/{x)A 



ik{x) 5n{x) 



U 



1 



k{x)v{x)h 



U 



ik{x) 6fi{x) 2 
1 5 k{x)v{x)h 



, .X{k{x),^{x),v{x),t) . (4.39) 

_ik[x) dfi{xj 2 J J J 

Now the comparison with the results of Q may be easily understood. Let us stick to quantum 
field theory for definiteness. The quantum characteristic functional which is a generating 
functional for correlation functions of the $ field coincides with the generating functional 
considered in ^ 

Z{n', u\ t) = Tr {p{t) exp | / (fx p!{x)(j){x) + u'{x)tt{x) 



(4.40) 

after rescaling the parameters and u to ^' = ikfj, and i'' = iku (of course, the same holds for 
quantum mechanics). Consequently, the evolution equations for the characteristic functional 



may be seen to be equal to those found in Q for the generating functional ( 4.40 ) provided the 
parameters /i' and u' are rescaled as specified. 

Going back to the initial remark of this section, we may conclude that, the quantum 
characteristic function and its evolution equation (or the generating functional in ( [4.40D ) are 
more interesting from an operative point of view as they determine the correlation functions 
and their time evolution. On the other hand the introduction of the MDF is both relevant and 
necessary from a theoretical point of view. In facts it allows a unified description of classical and 
quantum phenomena in terms of probability distributions obeying different evolution equations. 
Also it justifies the introduction of the X variable, as a variable taking values in an ensemble 
of reference frames (1.2), in view of the invertibility of the MDF in terms of the density matrix. 
This seems to us the profound motivation for introducing such a combination of phase space 
variables in the quantum characteristic functional and in the generating functional ( 4.40| ) . 
We stress once again that X and its field analogue ^ are, for each couple (^, z^), configuration 
space variables in the transformed reference frame labelled by (/x, v). 
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5 Conclusions 



In this letter we have presented an extension of the probabihstic representation of quantum 
mechanics to quantum field theory. In this framework classical and quantum phenomena, both 
statistically described, only differ by the evolution equations of the distributions of probabilities 
for the relevant observables. Quantum observables are described by a distribution of proba- 
bility, the MDF, and the time evolution by an integro-differential equation for the MDF. We 
recently addressed the problem of finding the Green's function for the time-evolution equation 



of the MDF |1C]. The problem was solved for quadratic Hamiltonians, and a characteriza- 
tion of such a propagator in terms of the time-dependent invariants of the system was found. 
This propagator represents the transition probability of the system from a quantum state to 
another. Thus, a generalization to quantum field theory would be interesting in our opinion, 
and is presently under consideration. Another promising application of the probabilistic point 
of view is suggested in where it is used to study the approach to equilibrium of non equi- 
librium quantum field theories. An extension to relativistic quantum field theory would be 
also interesting, though it poses problems of interpretation which are not understood at the 
moment. 
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